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Abstract 

Coherent dynamics of atomic matter waves in a standing-wave laser field is studied. In the 
dressed-state picture, wave packets of ballistic two-level atoms propagate simultaneously 
in two optical potentials. The probability to make a transition from one potential to 
another one is maximal when centroids of wave packets cross the held nodes and is given 
by a simple formula with the single exponent, the Landau-Zener parameter re. If re ^> 1, 
the motion is essentially adiabatic. If re <C 1, it is (almost) resonant and periodic. If 
re — 1, atom makes nonadiabatic transitions with a splitting of its wave packet at each 
node and strong complexification of the wave function as compared to the two other 
cases. This effect is referred as nonadiabatic quantum chaos. Proliferation of wave 
packets at n ~ 1 is shown to be connected closely with chaotic center-of-mass motion in 
the semiclassical theory of point-like atoms with positive values of the maximal Lyapunov 
exponent. The quantum-classical correspondence established is justified by the fact that 
the Landau-Zener parameter k specifies the regime of the semiclassical dynamical chaos 
in the map simulating chaotic center-of-mass motion. Manifestations of nonadiabatic 
quantum chaos are found in the behavior of the momentum and position probabilities. 

Keywords: cold atom, matter wave, quantum chaos 



1. Introduction 

The mechanical action of light upon neutral atoms has been comprehensively studied 
since the pioneer works of Lebedev, Gerlah and Stern, Kapitza and Dirac and Frisch. 
The light pressure force provides optical cooling and trapping of atoms [1]. In the last 
two decades, cold atoms in standing-wave optical helds have been used to study quantum 
chaos. The proposal Q to study atomic dynamics in a far-detuned modulated standing 
wave made atomic optics a testing ground for quantum chaos. A number of impressive 
experiments have been carried out in accordance with this proposal 0, 0, @| • New pos- 
sibilities are opened if one works near the atom-held resonance where the interaction 
between the internal and external atomic degrees of freedom is intense [f| 0, H| • 

Dynamical chaos in classical mechanics is a special kind of random-like motion with- 
out any noise and/or random parameters. It is characterized by exponential sensitivity 
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of trajectories in the phase space to small variations in initial conditions and/or control 
parameters. Such sensitivity does not exist in isolated quantum systems because their 
evolution is unitary, and there is no well-defined notion of a quantum trajectory. Thus, 
there is a fundamental problem of emergence of classical dynamical chaos from more 
profound quantum mechanics which is known as quantum chaos problem and the related 
problem of quantum-classical correspondence. In a more general context it is a problem 
of wave chaos. It is clear now that quantum chaos, microwave, optical, and acoustic chaos 
S fiol [ill [H| have much in common. The common practice is to construct an analogue 
for a given wave object in a semiclassical (ray) approximation and study its chaotic prop- 
erties (if any) by well-known methods of dynamical system theory. Then, it is necessary 
to solve the corresponding linear wave equation in order to find manifestations of classi- 
cal chaos in the wave-field evolution in the same range of the control parameters. If one 
succeeds in that the quantum-classical or the wave-ray correspondence are announced to 
be established. 

In atom optics [l3| one quantizes both the atomic internal and translational degrees 
of freedom. The atom is treated as a wave packet which undergoes deformations in the 
process of exchange of energy and momentum quanta with a light wave. Quantization 
of the translation motion provides an entanglement of the internal and external degrees 
of freedom. Any changes in the form of the wave packet will affect the internal state 



of the atom and vice versa [14|, [15|- The optical field provides a tool to manipulate 



the atomic matter waves. In atom optics the Schrodinger equation for the probability 
amplitudes constitutes a linear infinite-dimensional dynamical system which is governed 
by an external force if the field is treated as a classical wave. 

In the semiclassical approximation, atom with quantized internal dynamics is treated 
as a point-like particle with the Hamilton-Schrodinger equations of motion constituting 
a nonlinear dynamical system. A number of nonlinear Hamiltonian and dissipative dy- 
namical effects have been found with such a system including chaotic Rabi oscillations, 
chaotic atomic transport, dynamical fractals, synchronization, chaotic walking, and Levy 
flights [l|| E3, EM 0, Hfl 21 1 . Similar and new effects have been found numerically and 
described analytically with two- level atoms in a losseless cavity with a single quantized 



mode in the framework of the Jaynes-Cummings model [22j, |23|, l24|, [25| and the Tavis- 
Cummings model [26]. It has been shown that the coupled atom-field dynamics in a 
cavity can be unstable under appropriate conditions in the absence of any kind of in- 
teraction with environment. This kind of quantum instability manifests itself in fractal 
chaotic scattering of atoms |22l l23l [2J, |25[, in strong variations of reduced quantum 



purity and entropy 24 , 2^, 26|. correlating with the respective maximal Lyapunov expo- 
nent, and in exponential sensitivity of fidelity of quantum states to small variations of 
the detuning [m I25I] . 

The main aim of this paper is to establish a kind of the quantum-classical correspon- 
dence in transport properties of point-like atoms and atomic matter waves moving in a 
standing-waves field. Is the coherent evolution of the atomic matter waves really com- 
plicated in that range of the control parameters where the corresponding center-of-mass 
motion has been shown to be chaotic? 
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2. Wave-packet motion in a standing light wave 



The Hamiltonian of a two-level atom, moving along a one-dimensional classical standing 
wave laser field, can be written in the frame rotating with the laser frequency w/ as 
follows: 

P 2 h 

H = + -j{^a - u f )a z - htt (a- + a + ) cosk f X, (1) 

where a±, z are the Pauli operators for the internal atomic degrees of freedom, X and P 
are the atomic position and momentum operators, ui a and Q are the atomic transition 
and Rabi frequencies, respectively. We will work in the momentum representation and 
expand the state vector as follows: 

|*(*)) = J[a(P,t)\2) + b(P,t)\l)}\P)dP, (2) 

where a(P,t) and b(P,t) are the probability amplitudes to find atom at time t with the 
momentum P in the excited, |2), and ground, |1), states, respectively. After some algebra 
one gets the normalized Schrbdinger equation for the probability amplitudes [hi ] 

ia(p) = \{u r p 2 - A)a(p) - -[b(p + 1) + b(p - 1)], 

I I (3) 

ib(p) = -Kp 2 + A)b(p) - - [a(p + 1) + a(p - 1)], 

where the dot denotes differentiation with respect to dimensionless time r = fit, p = 
P/hkf, and x = kfX. The normalized recoil frequency ui r = hk 2 /m a fl and the atom- 
field detuning A = (ujf — u) a )/Q, are the control parameters. 

The probability to find an atom with the momentum p at the moment of time r is 
V(p,r) — \a(p, t)| 2 + \b(p 7 t)| 2 . The internal atomic state is described by the following 
real- valued combinations of the probability amplitudes: u q {r) = 2 Re f dp[a(p,r)b* (p, r)], 
v q (r) = — 21m J dp[a(p,r)b* (p, t)], z q (r) = J dp[\a(p, r)| 2 — \b(p, r)| 2 ], which are expected 
values of the synchronized (with the laser field) and a quadrature components of the 
atomic electric dipole moment (u q and v q , respectively) and the atomic population in- 
version, z q . Varying the value of the Rabi frequency f2, we can change the value of the 
dimensionless recoil frequency u r with the same atom. Working, say, with a cesium atom 
(m a = 133 a.u., \ f = 852.1 nm, and f rcc ~ 2 KHz), we get uj r = 10" 5 at fi = 100 MHz. 



We will interpret the wave-packet motion in the dressed-state basis [l3j, [27 1 

| + ) A = |2) sin 6 + |1) cos 6, | - ) A = |2) cos6 - |1) sin6, (4) 
where Q is the mixing angle 



tan9=-^ J(—^—\ +1. (5) 

2 cos x y \ 2 cos x J 

These states are eigenstates of an atom at rest in a laser field with the eigenvalues of the 
quasienergy 

E ( ^ ] = ±\J^- +cos 2 x. (6) 
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The probability amplitudes to find the atom at point x in those potentials are, respec- 
tively 

C+(x) = a(x) sin 9 + b(x) cos 9, C-(x) — a(x) cos 9 — b(x) sin 9, (7) 

where the amplitudes in the bare-state basis a(x) and b(x) may be computed in the 
position representation with the help of the Fourier transform 

/oo />oo 
dp'e zp ' x a(p'), b(x) = const / dp' e ip ' x b(p'). (8) 
-oo J —oo 

Let us assume that we are able to prepare an atom exactly in one of its dressed states, 
| + ) A or | — ) A . Then the atom will move in one of the potentials, E^ or \ along 
a single trajectory. In quantum mechanics, there is a nonzero probability to make a 
transition to another potential. To estimate this probability we write the Hamiltonian 
of the internal degree of freedom of a two- level atom in the basis | ± ) A 

flint = &z cos x + —a x . (9) 

Let us linearize the cosine in the vicinity of a node of the standing wave and estimate 
a small distance the atom makes when crossing the node as follows: Sx = uj r |p n ode | t 
[l3j | . The quantity w r |p no do| is a normalized Doppler shift for an atom moving with the 
momentum [p n ode|, i.e., lod = <^r\Pnode\ = kf\v no a e \/n. The Schrodinger equation for the 
probability amplitudes C± (x) in the position representation can be written in the form 
of the second-order equation 



C+(x) + 

The asymptotic solution of Eq. (fTTJl) 



A 2 

iuj D + — + (ujdt) 2 



C+(x) = 0. (10) 



Plz = exp(-/s), (11) 

gives the probability to make a nonadiabatic or Landau-Zener transition from one of the 
nonresonant potentials to another one specified by the Landau-Zener parameter 

A 2 

K = 7T . (12) 

(jJd 

There are three regimes of atomic motion. 

1. k> 1. The probability to make the transition is exponentially small even when an 
atom crosses a node. The evolution of the atomic wave packet is adiabatic in this 
case. 

2. k <C 1. The distance between the potentials at the nodes is small and the atom 
changes the potential each time when crossing any node with the probability close to 
unity. In the limit case A = 0, the atom moves in the resonant potentials. 

3. k ~ 1. The probability to change the potential or to remain in the same one, upon 
crossing a node, are of the same order. In this regime one may expect a proliferation 
of components of the atomic wave packet at the nodes and complexification of the 
wave function. 
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3. Simulation of ballistic wave-packet propagation 

We simulate the evolution of a Gaussian wave packet with the variance in the mo- 
mentum space, (Tp = 50, po = 10 3 , xq = and cu r = 10 -5 . The initial average kinetic 
energy, u! r p 2 /2 = 5, is greater than the depth of the potential wells, so the atom will 
move ballistically along the positive direction of the standing-wave axis. 

To study all the regimes of the wave-packet motion, we simulate Eqs. © at different 
values of the Landau-Zener parameter re (fT2]) . The normalized Doppler shift cup nearby a 
node of the standing wave is estimated to be ljd ~ w r po = 0.01. If we choose, say, A = 0.3 
we get the first case in our nomenclature, /t» 1, with exponentially small probability of 
nonadiabatic transitions. The wave packet, initially prepared in the ground state which 
is a superposition of the dressed states with approximately equal weights, splits from 
the beginning (Fig. [1^) into two components, | + ) A and | — ) A , each of which moves in 

its own nonresonant potential, or £" A ' . We really do not observe in Fig. [TJi any 

splitting at the nodes, and the motion of the wave packet at A = 0.3 is adiabatic and 
practically periodic. If k <C 1 (the motion near the resonance), one expects to observe 
the periodic motion in the two resonant potentials simultaneously without any splitting 
.14]. 

At A = 0.1, we get re ~ 1 and expect nonadiabatic transitions at the nodes of the 
standing wave in accordance with formula (1111) . The initial ground state |1) is now 
a superposition of the dressed states with practically the same weights. The initial 
bifurcation is accompanied by splittings (see Fig. [Tb) that can be proved to occur at the 
nodes of the standing wave. Let us start to analyze the wave-packet motion with the 
| + ) A -component (the upper curve in the figure). The first splitting occurs at r} + ' ~ 150. 
It is easy to prove that it is the moment of time when the centroid of the | + ) A -component 

crosses the first node at x = ir/2: t[ + ^ = 7r/2w,.pQ + L ' ) ~ 150, where Pq^ is the average 
momentum of the centroid of the | + ) A wave packet between x = and x = n/2. Thus, 
the wave packet, crossing the node, splits into two parts. The first one prolongs its 
motion in the potential E^ after passing the point x — n/2. It is the lower curve 
in Fig. QJ) starting at t± ~ 150. The corresponding packet slows down because this 
component loses its kinetic energy going up to the top of the potential E^ . As to the 
second trajectory (the upper curve starting at t[ ), it appears due to the nonadiabatic 
transition to the potential E^ '. That is why it accelerates from the beginning and 
reaches its maximal velocity at x — tt. The | — ) A -component (the lower curve starting 

at r = 0) splits at the first node at t[ ^ = it j2u r p\ 1 ~ 156. In course of time both the 

components split at every node of the standing wave at the moments r« that can be 
estimated with the simple formula 

"rpf-iA^ = ( 2n - l )\ > « = 2, 3, . . . , (13) 

where „ is an average momentum of the centroid of the corresponding component 
between the (n — l)-th and rt-th nodes. Such a proliferation at the nodes means a 
complexification of the atomic wave function both in the momentum and position spaces 
as compared to the adiabatic and resonant cases. 

Now we go to the position space and compute the probability \C(x, r)| 2 = |C_ (x, r)| 2 + 
\C + (x,t)\ 2 to be at point x at time r. In Fig. [5] we show the result of simulation in the 
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case of adiabatic and nonadiabatic motion at A = 0.3 and A = 0.1 corresponding, respec- 
tively, to Fig.[T^, and b in the momentum space. It is a plot of the position probability in 
the frame moving with the initial atomic velocity uJ r Po — 0.01 where the slope straight 
lines mark positions of the nodes of the standing wave in the moving frame. At A = 0.3, 
the evolution is simple without any transitions at the nodes (Fig. [5^) . The splitting of 
the total probability \C(x, r)| 2 is caused by the initial bifurcation of the wave packet due 
to its bipotential motion. 

The situation is cardinally different when we work in the regime with nonadiabatic 
transitions at the held nodes (A = 0.1). Splitting at the nodes in the momentum space 
(see Fig[T}D) manifest itself in the position space in Fig. ^jp. In this case one observes 
visible changes in the proba6bility |C(a;)| 2 exactly at the node lines. It is a clear evidence 
of the nonadiabatic transitions that occur in the specific range of the control parameters, 
k ~ 1. This results in a proliferation of components of the wave packet at the nodes and, 
therefore, a complexification of the wave function both in the momentum and position 
spaces. 



4. Quantum-classical correspondence and nonadiabatic quantum chaos 

In this section we compare the quantum results, obtained in the preceding sections, 
with those obtained for the same problem but in the semiclassical approximation when 
the translational motion has been treated as a classical one @,0> B 12, 3- Coherent semi- 



classical evolution of a point-like two-level atom is governed by the Hamilton-Schrodinger 
equations with the same normalization as in the quantum case 

x = ijj r ri. v = — it sin ,-r, u — Av, 

(14) 

v — — Au + 2zcos:r, z — — 2i>cosa;, 

where 

u = 2Re(a &5), v = -2lm{a b*), z = \a \ 2 ~ \b a \ 2 (15) 

are the atomic-dipole components (u and v) and population-inversion (z), and clq and bo 
are the complex- valued probability amplitudes to find the atom in the excited and ground 
states, respectively. The system (fl4|l has two integrals of motion, the total energy 

CJ r 2 A 

H = —p — ucosx — —z, (16) 

and the length of the Bloch vector, u 2 + v 2 + z 2 = 1. 

Equations Q14p constitute a nonlinear Hamiltonian autonomous system with two and 
half degrees of freedom and two integrals of motion. It has been shown in Ref. [19j to 
have positive values of the maximal Lyapunov exponent A in a wide range of values of the 
control parameters and initial states. This fact implies dynamical chaos in the sense of 
exponential sensitivity to small changes in initial conditions and/or control parameters. 
The result of computation of the maximal Lyapunov exponent in dependence on the 
detuning A and the initial Doppler shift cjd = ui r po is shown in Fig. [3] at uj r = 10~ 5 . In 
white regions of the plot the values of A are almost zero, and the atomic motion is regular 
in the corresponding ranges of A and u>d ■ In shadowed regions positive values of A imply 
unstable motion. At exact resonance, we get A = because at A = the semiclassical 
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equations of motion (fl4| become integrable due to an additional integral of motion, 
u = const. We stress that the local instability produces chaotic center-of-mass motion 
in a rigid standing wave without any modulation of its parameters in difference from 
the situation with atoms in a periodically kicked optical lattice [!, 0, [Bj] . In dependence 
on the initial conditions and the parameter values, an atom may oscillate in a well of 
the lattice or it may have enough kinetic energy to overcome the potential barrier. In 
some cases the center-of-mass motion resembles a random walking. It means that an 
atom in a deterministic standing-wave field alternates between flying through the lattice, 
and being trapped in its wells. Moreover, it may change the direction of motion in a 
random-like way (see Ref. 17, ljl] for coherent Hamiltonian dynamics and Refs. 2(J 21 1 



for a dissipative one with spontaneous emission included). 

It follows from (fT4|) that the translational motion is described by the equation for a 
nonlinear physical pendulum with the frequency modulation 

x + uj r u(r) sin a; = 0, (17) 

where u is a function of all the other dynamical variables. It has been shown in Ref. [l9[ 
that the regime of the center-of-mass motion is specified by the character of oscillations of 
the component u of the Bloch vector. In a chaotic regime sudden "jumps" of the variable u 
occur when an atom crosses the field nodes. Figure^ demonstrates more or less periodic 
oscillations of u at the detuning value A = 0.3 at which the corresponding maximal 
Lyapunov exponent is zero (Fig. [3]). In the chaotic regime at A = 0.1 u demonstrates 
shallow oscillations interrupted by jumps of different amplitudes upon crossing the nodes 
(Fig.Sb). 

Approximating the variable u between the nodes by constant values, the following 
stochastic map has been constructed in Ref. [l9| 

u m — sin(0sin!y9 m + arcsinit m _i), (18) 



where G = y 7rA 2 /w r p n odc is an angular amplitude of the jump, u m is a value of u just 
after the m-th node crossing, ip m are random phases to be chosen in the range [0, 2ir], 
and Pnodc = \/2H/uj r is the value of the atomic momentum at the instant when the 
atom crosses a node (which is the same with a given value of the energy H for all the 
nodes). With given values of A, uj r and p n ode, the map (|18|) has been shown numerically 
to give a satisfactory probabilistic distribution of magnitudes of changes in the variable 
u just after crossing the nodes. The stochastic map (TT51) is valid under the assumptions 
of small detunings (|A| <C 1) and comparatively slow atoms (|w r p| <C 1). Furthermore, 
it is valid only for those ranges of the control parameters and initial conditions where 
the motion of the basic system JTJJ is unstable. For example, in those ranges where all 
the Lyapunov exponents are zero, u becomes a quasi-periodic function and cannot be 
approximated by the map (see Fig. 0^) . 

The key result in the context of the quantum-classical correspondence is that the 
squared angular amplitude of the map (|18|) is exactly the Landau-Zener parameter (|12p , 
i.e., 2 = k. Rewriting the map (|18[) for arcsinu m , one gets 

arcsinu m = y/Kshitp m + arcsinit m _i, (19) 



where the jump magnitude does not depend on a current value of the variable. The map 
(fT§|) visually looks as a random motion of the point along a circle of unit radius (see Fig. 4 
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in Ref. [19[). If K ~ 1, then the internal atomic variable arcsinu m just after crossing 
the m-th node may take with the same probability practically any value from the range 
[— 7r/2, 7r/2]. It means semiclassically that the momentum of a ballistic atom changes 
chaotically upon crossing the field nodes. In accordance with the quantum formula (|12[) . 
the corresponding atomic wave packet makes nonadiabatic transitions when crossing the 
nodes and splits at each node (see Figs. QJ> and^)). As the result, the wave packet of 
a single atom becomes so complex that it may be called a chaotic one in the sense that 
it is much more complicated than the wave packets propagating adiabatically. Thus, 
nonadiabatic wave chaos and semiclassical dynamical chaos occur in the same range of 
the control parameters and are specified by the same Landau-Zener parameter n ~ 1. In 
two limit cases with k <C 1 and k> 1 both the semiclassical and quantized translational 
ballistic motion are regular. 

In quantum mechanics there is no well-defined notion of a trajectory in the phase 
space and, hence, the Lyapunov exponents can not be computed. In quantum mechanics 
there is no exponential sensitivity to small variations in initial conditions because the time 
evolution of an isolated quantum system is unitary, and the overlap of any two different 
quantum state vectors is a constant in course of time. Moreover, quantum phase space 
is discrete due to the Heisenberg uncertainty principle unlike continuous classical phase 
space. Namely the continuity of the classical phase space provides a possibility of chaotic 
mixing which exploits more and more fine structures in the classical phase space in course 
of time whereas the quantum evolution stops to do that over a rather short Ehrenfest 
time. 

The semiclassical (fl"4)) and quantum <j3j) equations of motion look very different. The 
semiclassical ones constitute a five-dimensional nonlinear dynamical system of ODEs 
with two integrals of motion that has been shown to be chaotic in a certain range of 
control parameters with exponential sensitivity to small variations to initial conditions 
fl9^ . The quantum ones constitute an infinite-dimensional set of linear equations. It is 
not evident a priori that their solutions might demonstrate a kind of correspondence 
in the same range of the control parameters. Nevertheless, a sort of quantum-classical 
correspondence both in regular and chaotic regimes of the center-of-mass motion has 
been found. It should be stressed that this correspondence manifests itself in behavior 
of the quantum Bloch variable u in the semiclassical equations of motion (1141) . 

However, the quantum-classical correspondence is not and could not be absolute 
because the Planck constant is equal to 1 with our normalization. It cannot tend to be 
zero in order to achieve a classical limit as it could be done with an effective Planck 
constant (see, for example, Refs. @, depending on the system's parameters. We work 
in this sense in a deep quantum regime. The quantum-classical dualism with cold atoms 
resembles the wave-ray one in a classical wave motion. In the context of this paper we 
might compare ray-like trajectories of atoms with their wave-like motion. 

To illustrate correspondence and difference that inavoidably appears when comparing 
quantum evolution with the classical one (that is only an approximation to the quantum 
one), we compute with Eqs. ([14")) the evolution of a Gaussian distribution over classical 
momentum p and position x with the same parameter's values as in simulation of the 
wave-packet propagation shown in Figs. [T] and [2] In accordance with the Lyapunov 
map in Fig. [3j one expects a regular center-of-mass motion at the detuning A = 0.3 and 
a weakly chaotic one at A = 0.1. In Figs. [SJi and b evolution of classical momenta is 
shown for the regular (A = 0.3) and chaotic (A = 0.1) regimes of the center-of-mass 
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motion. Visible spreading in p with chaotically moving atoms, as compared to regularly 
moving ones, is one of the signs of classical dynamical chaos. Figure [1] demonstrates 
similar spreading of the momentum probability distribution of a Gaussian wave packet 
with nonadiabatic transitions at A = 0.1 (Fig. [1}d) as compared to the adiabatically 
moving wave packet at A = 0.3 (Fig. QJ,). The difference between the classical and 
quantum evolution is also evident: the semiclassical equations of motion (TH|) are not 
able to simulate the splitting of wave packets due to purely quantum effect of motion in 
two optical potentials simultaneously 

In Figs. and b we plot, respectively, regular and chaotic trajectories in the frame of 
reference moving with the initial atomic velocity uJ r Po — 0.01. The bundle of chaotically 
moving atoms in Fig. [SJd diverges in a short time significantly as compared to the regular 
one in Fig. [6^. This property can be used to detect chaotic scattering in a real experiment 



with atoms crossing a standing laser wave [151 ] . As to quantum motion in the position 



space, it is evident that the wave packet with nonadiabatic transitions (Fig. ^p) becomes 
much broader in course of time due to splitting at the standing-wave nodes than the 
adiabatic wave packet in Fig. [3^ resembling a broadening of the bundle of chaotic point- 
like atoms (Fig. |6}d) as compared to the regular one (Fig. However, we do not 
observe any splitting of the classical bundles because a classical trajectory simulates only 
the motion of the centroid of a quantum wave packet and cannot simulate of course its 
splitting due to purely quantum effect of motion in two optical potentials simultaneously. 
There is only one optical potential in the semiclassical approximation. 



5. Conclusion 

We have studied coherent dynamics of ballistic atomic wave packets in a one-dimensional 
standing-wave laser field. The problem has been considered in the momentum represen- 
tation and in the dressed-state basis where the motion of a two-level atom was interpreted 
as a motion in two optical potentials. The character of that motion has been shown to 
depend strongly on the value of the Landau-Zener parameter k (|12p . If k ^ 1, then the 
probability of transitions from one of the potential to another one, which is described by 
the Landau-Zener formula dill) , is exponentially small. Under such a condition, atoms 
move in the adiabatic regime. If k <C 1, the formula (|llj) gives almost unity probability 
to change the potential when crossing the nodes. In the intermediate 1, the 

probabilities for an atom to change or not to change the nonresonant potential, when 
crossing a node, are of the same order. The corresponding nonadiabatic transitions man- 
ifest themselves as a splitting of the atomic wave packets in the momentum space when 
their centroids cross the nodes. This nonadiabatic quantum chaos occurs exactly in the 
same range of the detuning and the Doppler shift where the semiclassical dynamics has 
been shown to be chaotic. It is remarkable that the same Landau-Zener parameter k 
specifies both semiclassical and quantum chaos with ballistic atoms in a deterministic 
optical lattice. 

We hope that the results obtained can be used to study manifestations of quantum 
chaos with Bose-Einstein condensates in optical lattices [28( with coupled degrees of 
freedom. From the theoretical point of view, the dynamics of condensates of ultracold 
atoms is described correctly by the Gross-Pitaevskii equation which is a kind of a non- 
linear Schrodinger equation with possible chaotic solutions. Experimentally, one of the 
possibilities is to prepare two Bose-Einstein condensates in different internal states [29| . 
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Another possibility can be realized with a Bose-Einstein condensate in an optical lattice 
subject to a static tilted force 3^, 31 , 32|. Viewing transitions between the Bloch bands 
of a condensate in such a tilted optical lattice as a two-state problem, we get a meso- 
scopic quantum system with coupled different degrees of freedom (I am thankful to an 
anonymous referee for that comment). 
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Figure 1: (Color online) Momentum probability distribution V(p, t) of a Gaussian wave packet vs time 
with po = 1000, <Tp = 50, and lo t = 10~ 5 at (a) A = 0.3, adiabatic motion, and (b) A = 0.1, motion 
with nonadiabatic transitions. The color codes the values of P(p, r). 
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Figure 2: (Color online) The position probability C(ai)j 2 in the moving frame of reference with the slope 
straight lines marking positions of the nodes, (a) Adiabatic motion in the position space at A = 0.3. 
(b) Wave-packet propagation in the position space with nonadiabatic transitions at the field nodes. 




Figure 3: Maximal Lyapunov exponent A vs atom-field detuning A and the initial Doppler shift uj£> 
oJrPo- Color codes the values of A. 
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Figure 4: Semiclassical evolution of the atomic-dipole component u in (a) regular (A = 0.3) and (b) 
chaotic (A = 0.1) regimes of the ballistic motion of a point-like atom. 
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Figure 5: Atomic trajectories in the momentum space computed with the classical Gaussian distribution 
at the same parameter's values as in simulation with Gaussian wave packets, (a) Regular center-of-mass 
motion (A = 0.3) corresponding to adiabatic quantum motion in Fig.^i and (b) weakly chaotic motion 
(A = 0.1) corresponding to nonadiabatic quantum motion in Fig. [lp. 




Figure 6: Classical atomic trajectories in the moving frame of reference, (a) Regular bundle (A = 
0.3) corresponding to adiabatic quantum motion in Fig. [2^, and (b) weakly chaotic bundle (A = 0.1) 
corresponding to nonadiabatic quantum motion in Fig. [2p. 
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